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Abstract 

We consider the generalized exclusion statistics in the Kondo prob- 
lem. The thermodynamic Bethe ansatz equations have been used for a 
multicomponent system of particles obeying the generalized exclusion 
principle. We have found a relation between the derivative of the phase 
shift of the scattering matrix for Fermi particles and for particles charac- 
terized by generalized exclusion statistics. We show that the statistical 
matrix in the Kondo problem has a universal form in high and low tem- 
perature limits. 

PACS numbers: 05.30.-d, 72.15.Qm, 11.25.Hf 

1. Introduction 

Generalized exclusion statistics (GES) has attracted much interest in recent papers [2 - 
8] A typical example of the application of generalized exclusion statistics is the Cologero- 
Sutherland model §, |l^, D. Bernard and Y.-S. Wu |l| studied the GES in the 
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thermodynamic Bethe ansatz (TBA) using equivalence of the Bethe ansatz (BA) equations 
and the ones obtained by the Haldane principle jl^. M. Wadati has shown ||l3| that a 
change in the statistics is determined by an appropriate choice of the phase shift of the 
scattering matrix. The generalization of exclusion statistics for multicomponent systems 
was provide in the paper js). 

The goal of this paper is to consider the exclusion statistics of excitations in the Kondo 
problem. To solve this problem we have applied the TBA for a multicomponent system of 
particles obeying the GES. We have found the universal behavior of the statistical matrix 
in high- and low-temperature limits. In these limits the statistical matrix is proportional 
to the Cartan matrix for the A„-algebra. 

The paper is organized as follows. In the second section we provide the main equations 
of the exclusion statistics theory to complete the discussion. In this section we have also 
derived the relation between the derivative of phase shift (DPS) of the scattering matrix 
and the statistical matrix in the framework of the TBA equations for the multicomponent 
system of particles. The third section deals with the application of this result to the Kondo 
problem. We show that distribution function and statistical matrix in this problem have 
the universal form in high and low temperature limits. 



2. Exclusion Statistics Equations 

A change in the number D^i of the vacant states due to the addition of the number Nf}j 
of the particles, according to Haldane |^ is defined as 

(1) 

Here gai,pj is the matrix of statistical interaction. The indices a {a = 1,2, ...M) and 
i correspond to the internal and the dynamical degrees of freedom, respectively. The 
solution of Eq.(|l|) has the form: 

i^a. = -^5az,«iV^,+-Da., (2) 

where D^,^ is the number of vacant states of the ai-th type without particles. The number 
of holes Dai determines the statistical weights W as follows 
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In the specific cases of gai,f3j = and gai,i3j = ^ap^ij Eq.(^) yields the well-known expres- 
sion for the statistical weights of Bose and Fermi particles. 

The distribution function Uai is defined usually ^ ^ in the following way Uai = 
Nai/C^i- This is not convenient for systems with internal degrees of freedom. For example, 
in the hierarchical basis of the states in the fractional quantum Hall effect Uai = oo if 
a = 2, 3, ... because = for spin degrees of freedom. The definition of the distribution 
function in the form 

riai = -TT- (4) 

with 

Gc.^ = 0°, +No,^-Y, 9c.^,0JNp, (5) 

is more convenient because it takes into account the influence of the statistical interaction 
on the number of the states. The equilibrium distribution function be found in 

this case from the extremum of the grand partition function as a solution of the following 
equations: 

— n [1 ~ «/3jf = exp{(e,, - Ai„)/T} . (6) 

Here eai is the energy and /Iq is the chemical potential for the particles of type a. The 
distribution function Uai determines the free energy 
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(7) 



as well as the value of the entropy S = InW. 

The interection between quasiparticles from the TBA point of view is encoded in the 
phases Qai,pj of the dynamical scattering matrix Sai^pj — — exp^—iQai.pj)- The statistical 
properties expressed by the statistical matrix gai,pj depend on the DPS of the scattering 
matrix || 

Let us consider the TBA equations for the set of multicomponent particles obeying the 
GES. Quantizing a gas of such particles on a circle of length L requires that the momentum 
kai of the ai-th particle satisfies the following condition: 

N 

exp{ikML} J2 So.AQ^ ~ 0,) = 1 . (8) 

13] 
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The momentum and the energy of the particles are parametrized by the rapidity 9. Going 
to the hmits L ~> oo and N oo with the finite value of N/L and taking the derivative 
of the log of Eq.(|) yield 



dka{9) 



M 

E 

(3=1 



where 



(9) 



(10) 



is DPS. The function pp{9) in Eq.(^) is the density of the particles of /3-type, Qaid) is the 
density of the states. 

The information about the statistical properties of the system is contained in the 
distribution functions and in the entropy S. In the framework of the thermodynamic 
Bethe ansatz JlTt , we have neither information about the ground state energy, nor about 
the structure of low- lying excitations with the energy e" (6'). They are the objects of the 
traditional Bethe ansatz approach. The equilibrium state at temperature T is obtained 
by minimizing the free energy F — E — TS, where the energy of the system is 



M 



(11) 



The variation of F with respect to pa yields the following equations for the dressed 
energy €ai9): 



T 
2^ 



M 
[3=1 



1 



l-np{e') 



de' . 



(12) 



Here the function Ua relates to the functions pa and qa in Eq.(^ as follows: = Palla- 
Let us assume that the particles are fermions, i.e., 5/3j,ai = SafsSij. We see from Eq.(^ 
that 

1 , ^ 

n„ = p . (13) 



1 + exp 



After substituting this expression into Eq.(|l^), we have the standard TBA equations 



M 



13=1 



1 + exp 



T 



de' (14) 
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for fermions (superscript "/" denotes the Fermi statistics). 

Each statistics corresponds to the specific value of DPS $q^(0) in Eq.(p^. The tran- 
sition to the Fermi statistics leads to the new value ${^(6') of this function. Assuming 
that the function n^(6') in Eq.jT^) coincides with Eq.(p^ we can find from Eqs.jT^), (|lj) 
the relation between DPS ^ai3{(^) for the arbitrary statistics and for the Fermi statistics. 
This relation has the form 



$(61 - 6') ^^f{e-0')~ 2TTd{0 - 6') [I-G]. 
where X is the unit matrix and 

\ 

.921 .922 .923 

g 



931 932 933 



V 



(15) 



(16) 
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is the statistical matrix. It is easy to see from Eq.(p^ that for models with the function 
${q,(0 — 6') ^ S{9 — 6') one can find such statistical matrix Q which gives a zero value of 
the r.h.s. in this equation. In these models the phase of the scattering matrix has the 
structure of the step function. From the TEA point of view these systems look like a 
gas of non- interacting particles having the GES. This case is known ^ as the ideal 
exclusion statistics. In these models the correlations between particles can be transformed 
to the statistical interaction. The distribution function of excitations in the systems with 
the ideal exclusion statistics can be obtained from Eq.(|^) where the statistical matrix 
is now ga,i3 = Sa,p — 27*^Q/3 ^^'^ dressed energy of excitations coincides with bare 
energy e°. Note that the structure of Eq.(p|) looks like that for Eq.(p^ (after some simple 



transformation) because the function Ua satisfies Eq.(|l3| 



2. Generalized exclusion statistics 
Kondo problem. 



in 



the 



Let us consider the Kondo problem from the GES point of view. Bethe-ansatz equations 
for the isotropic s — d exchange model [n8| in the Kondo problem are 



exp{ikjL) = exp 



- i/2 



(17) 
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X^ + i/2Y fX-f + l/g + iS 



Xj—i/2j \\j + l/g — iSj — A/3 



n 

/3 



— Xf3 + i 



(18) 



These equations solve the problem of diagonalization of the s — d exchange hamiltonian 
with the impurity spin S and with the coupling constant I. Here the total number of (up- 
spin) electrons is denoted as N{P). The general solutions of Eqs.(|l7|), (|8|) have the form 
of n-strings according to the string hypothesis JT^ . The n-string is a set of n solutions 
given by 

A("'^"' ^X; + ^ - , J = 1, n . (19) 

Here A" is the real number and n is the order of the string. The distribution of the n-type 

particles (holes) density in the thermodynamic limit is Pn(A) {pl!^\x)). Assuming that 
the particles obey the Fermi statistics, the TEA equations for the function e„(A) have the 
form 



e„(A) = 



2eF , _i / X \ r T 

tan exp(7rA)d„i + 7^- / 

Ti" 2.7: J_, 



1 



cosh(A - A') 



In 1 



dX' . 



(20) 



Here and below the dressed energy is e„(A) = nn[p^''^(A)/p„(A)]. 

The solutions of these nonlinear integral equations describe the thermodynamic prop- 
erties of the s — ff-model. The external magnetic field H of the problem enters in the 
boundary condition as follows: 



(21) 



Note that this boundary condition means the condition of compensation of the internal 
magnetic field h = lim ( e„ / n) by external magnetic field H. This situation takes 

n — »oo 

place as well for some phase states in (2 + 1) _D systems. 
The spin free energy 



psp ^ _J^J. 



i— ^In (l + e'^(^)/^) dX 
h(7rA V / 



2 cosh(7rA) 



2 cosh(7rA + 1 /g 



ln(l + e^-(^)/^) dX 
•A + 1/q V J 



(22) 
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is expressed by the functions e„(A). The first term in the r.h.s. of Eq.(p^ corresponds 
to the spin free energy in the absence of impurity. The second term is the impurity 
contribution to the free energy. We will focus on the universal properties of the solutions 
of Eqs.® in the limits T ^ oo and T ^ 0. 

Comparing BA Eqs.(|20|) and TBA Eqs.([l^) one can see that the relation between them 
exists if 

eiiX) = -e,,{X) . (23) 

In other words, the energy of the particles in TBA Eqs.(^^ corresponds to the energy of 
the holes in BA Eqs.(pO|) of the Kondo problem. The a-type particle in the TBA equations 
is made corresponding to the n-th string solution of Eqs.(pO|). By changing the index n 
by a we can rewrite Eqs.(|2^) in the new notations as follows 



<^LW = —tan ^ exp(7rA)(5Qi - / 
TT 27r 



1 



cosh(A - A') 



In I 



dX' . 



(24) 



The boundary conditions have now the form = oo and lim e^/a = 

a — *oo 

and Eqs.(|lJ) we have 



and 



where 



$^(A- A') = 



2eF 



cosh(A - A') 



= tan ^ exp(7rA)(5ai , 



fO 1 
1 







1 0/ 



-H. FromEqs.(||) 
(25) 

(26) 
(27) 



Unlike in the TBA Eqs.(|T^, we know the function e° (A) in Eqs.(26) because it was found 
by the BA method. 

Let us assume now that the particles have the GES with any statistical matrix Q. In 
this case the matrix $ in Eqs.(p^) for the Kondo problem is 



$(A - A') 



cosh(A - A') 



-C-2TT6{X-X')[I-g] 



(28) 
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It is easy to see from the last expression that in the general case we cannot find a matrix Q 
to set the matrix $ to zero. However, it is possible to consider the particular cases of high- 
and low-temperature limits. The asymptotic behavior of the spin free energy determines 
the different values of the rapidity range required for the consideration of the high- and 
low-temperature limits. When T — > oo (weak coupling limit), the impurity free energy is 
given by 



imp = - — In I 1 + exp I - 



T 



(29) 



When T ^ (strong coupling limit). 



T 



Fimp = - -;t In 1 + exp 



^2S 



i+oo) \ 



T I 



(30) 



Therefore, in these temperature limits we have to consider ('^) l-^l ^ Taking the 
limit |A| ^ 00 in Eqs.(p4[) we have the equations for e^(±oo): 



T 



1 + exp 



eLi(±oo)^ 



T 



1 + exp - 



T 



(31) 



The boundary conditions are: lim e^(±(X))/a — —H, el{—oo) — 00 and e{(-|-oo) = 00. 

a — >c >Q 

The solution of these equations is [llSl 



-Tin 



/ sinh [^{a + 1) 
I sinh(^) 



(32) 



This result will be the same if we assume that the matrix $-^(A) in Eqs.(^J) is propor- 
tional to the 5- function, i.e., 



lim $^(A-A') = TrS{X-X')C. 



(33) 



It is clear because the function lim l/[cosh(A — A')] ~ 2/[exp(oo — A')] acts as the 5- 

X — ^00 

function at A' = 00. Therefore, in the GES case the matrix <I>(A) has the following form 



$(A - A') = 7r5(A - A') [C - 21 + 20] 



(34) 
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for these temperature limits. From the condition $ = of the ideal statistics we have the 
following form of the statistical matrix: 



g = 



vo 



-1 



0\ 

-1 

2 / 



(35) 



We see that it is proportional to the Cartan matrix for algebra, An- 

To find the distribution function Ha we need to know the structure of the low-lying 
excitations. It follows from Eq.(p4|) that ea{X) = for all a 7^ when A = —00 and for 
all a ^ I when A = +00. The structure of the equations for the distribution function of 
particles Ua, which can be obtained from Eq.(j^) with the statistical matrix (|35|), resembles 
that of the equations for eQ(|A| — 00) (^: 



1 



-1 = (1 - n^+i) (1 - 



The solution of Eqs.(|36|) (with the boundary condition lim ( — 
coincides with the solutions (p^) : 



- 1 



/T 



sinh(^(Q! 




at T ^ and 



sinh(^) 



sinh(^a) 
sinh(^) 



(36) 



exp(-aff/r)) 



(37) 



(38) 



at T — ^ 00. The spin free energy does not vary when we change the particle statistics. 

The function does not change in the TBA equations for multichannel Kondo effect 
p6| , p^ , |20t . But the function e5^(A) is changed. It depends on the number of channels, k, 
as 

= Sak exp ttA . (39) 
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This leads to the new solutions for e(, {X) in the low-temperature limit ( A 



+ 00): 



e{(+oo) 



Tin 
Tin ( 



sin(7r(a+l)/(fc+2) 
sin(7r/(fe+2)) 



sinh(^(Q + l-fe) 
sinh(^) 




a > k 



a < k 



(40) 



The solutions for the higher-temperature limit (A — > —00) are not changed. To find the 
solution ( ^o|) using the statistical matrix (|3^) we have to impose the additional boundary 
condition = 1 at a = fc. 



Let us discuss in conclusion what new insights into the thermodynamics of integrable 
systems (the Kondo problem in particular) are gained by considering these systems in 
terms of the GES. 

Using the GES principle we introduce an additional "parameter", that is, particle 
statistics determined by the form of the statistical matrix. If we suppose that the statistical 
matrix is arbitrary, we can write the TBA equations for a system of particles with any 
statistics. Each form of the statistical matrix in this case has a corresponding distribution- 
and the DPS-function. Therefore, we can write the TBA equations in a more convenient 
form. A successful choice of the statistical matrix may lead to a simpler form of the 
coupled equations for the dressed energy and for the distribution function. In particular, 
in the case of ideal statistics one can find a statistical matrix such that the solution of the 
TBA equations for the dressed energy coincides with the bare energies, and the equations 
for the distribution function repeat (up to transformation) the TBA equations for Fermi 
particles. 

Let us focus on the features of the Kondo problem which are studied using the GES 
principle. The consideration of spin excitations as quasiparticles in the Kondo problem 
leads to the conclusion that spin excitations correspond to the holes in the TBA approach. 
To find the distribution function in the system of particles with internal degrees of free- 
dom we need additional information besides the statistical matrix. We should know the 
detailed structures of the ground state j|] . The energy of the ground state shows itself 
as the boundary conditions for the equations determining the distribution function. The 
boundary conditions include the external magnetic field which compensates the internal 
"conformal" magnetic field of the system. 



Conclusion 
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The solution for the distribution function in high and low-temperature regions has the 
universal form. It is determined by the (/-deformed dimension [a + 1 ] , of irreducible 
representations of the quantum group ?7q (5^2) • Here [ x ] q = {q^ — q~^)/{q — q^^)- In the 
case of the single channel Kondo problem, the deformation parameter q = exp [H / 2T) 
is real. It depends on the external magnetic field and the temperature. In the multichannel 
Kondo problem (at a < fc), the deformation parameter q = exp[j7r/(fc + 2)] being the root 
of unity, is determined by the number, k of the channels. 

In summary, we used the TBA equations for a system of particles with internal degrees 
of freedom to consider the features of the GES in the Kondo problem. It is shown that 
the statistical matrix and the distribution function in the Kondo problem have a universal 
form in high- and low-temperature limits. 
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